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One of the most efficient approaches to control a multiple degree-of-freedom robot 
manipulator is the virtual decomposition control (VDC). However, the use of the re- 
gressor technique in the conventionnal VDC to estimate the unknown and uncertaities 
parameters present some limitations. In this paper, a new control strategy of n-DoF 
robot manipulator, refering to reorganizing the equation of the VDC using the time 
delay estimation (TDE) have been investigated. In the proposed controller, the VDC 
equations are rearranged using the TDE for unknown dynamic estimations. Hence, 
the decoupling dynamic model for the manipulator is established. The stability of the 
overall system is proved based on Lyapunov theory. The effectiveness of the proposed 
controller is proved via case study performed on 7-DoF robot manipulator and com- 
pared to the conventionnal Regressor-based VDC according to some evalution criteria. 
The results carry out the validity and efficiency of the proposed time delay estimation- 
based virtual decomposition controller (TD-VDC) approach. 
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1. INTRODUCTION 

In many robotic applications, the principal technical challenges arise in control implementations, spe- 

cialy when its subject to hight number of degree-of-freedom (DoF). Indded, the robotic systems with high DoF, 
can be modeled by a set of coupled highly nonlinear differential equations, with various uncertainties and dis- 
turbances, which increases the complexity of their control. 
A wide range of approaches have been proposed in the literature to control robot systems, ranging from lin- 
ear to nonlinear techniques, such as computed torque control (CTC), robust control, passivity based control, 
Lyapunov stability based robust control, sliding mode control (SMC) [1]-[7], Feedback Linearization , Back- 
stepping [8]-[12]. All these techniques are based on the traditional Lagrange-Euler formulation, which present 
inherent inconvenient, in complexity, then in computational burden [13]-[18]. 

More we have a high number of DoF in the system, more than the complexity of the dynamic model 
and the computation burden will increase [19], [20]; It was proportional to the fourth power of the DoF of the 
robot [21]. This problem limits the use of these algorithms and reduces the feasibility of the control system. 

Hence, to cope with the aforementioned problem, a novel theory based on virtual decomposition con- 
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trol (VDC) was proposed by [21] to solve the modeling and control problems of a multi-DoF robotic system. 
It is defined as adaptatif control approach [22]-[25]. In the VDC approach the entire system is decomposed 
virtually into subsystems ( single joint and single link). The dynamic interaction between two adjacent subsys- 
tems is handled by virtual power flow (VPF), which leads to prove the virtual stability. For each subsystem, a 
subcontroller is designed independently, while the stability of the global system is rigorously maintained [22], 
[24], [26]-[30]. The obtained dynamic equation of each subsystem is relatively simple, which can reduce the 
computation burden of the system. Compared to the Lagrangian dynamic model, the computation of VDC 
method is proportional only to the number of subsystems (DoFs). In order to define the full dynamics of the 
investigated robot, the VDC approach is based on regression technique. Indeed, the dynamic equation of each 
subsystem can be linearly parameterized in terms of a regressor matrix, and a unknown parameters vectors, 
using the required linear/angular velocity vector and its time-derivative. However, the regressor-based VDC 
control presents difficulties in practical implementation due to the complexity of the regressor matrix, espe- 
cially to estimate the uncertainties parameters [31], [32]. Indeed, the regressor matrix is known to be highly 
nonlinear, and its derivation is not unique and remain tedious, although the process is standardized which 
increase the computation complexity. 

In order to avoid the use of the regressor techniques, some alternative tools have been proposed in the 
literature [33]-[36], using the approximation-based adaptive control based on function approximation technique 
(FAT). In order to approximate the uncertainties parameter vectors of the dynamic model, the FAT technique is 
based on linear parameterization in the form of the weighting matrix and the (orthogonal) basis function matrix 
of the target matrix/vector. However, it presents some limitations. Indeed, the FAT is selected arbitrary without 
specific criteria, thus an approximation error can be produced. Furthermore, the estimation of the initial values 
of the weighting matrix is complex. Using the FAT for high DoF can be generated with a computational com- 
plexity of weighing and orthogonal matrices [33]. In order to propose a suitable solution to the aforementioned 
problems, this paper introduces a new nonlinear adaptive control strategy including the virtual decomposition 
control (VDC) [21] and time delay estimation (TDE) [37], known with its several advantages. The basic idea 
of the proposed time delay estimation-based virtual decomposition controller (TD-VDC) can be summarized 
as: 

a. The VDC is used as an efficient tool to handle the full-dynamics-based control problem of n-DoF robot ma- 
nipulator. This approach considers the dynamics of subsystems (rigid bodies and joints) to carry out a tracking 
trajectory, while guaranteeing the stability and convergence of the entire robotic system. 

b. The TDE is used to estimate simply and effectively, the unknown parameters vectors including uncertainties 
and external disturbances. It require a use of time-delayed information of the control torque inputs and state 
derivatives for each subsystems. 

c. Based on the aforementioned advantages of VDC and TDE approachs, the TDVDC is used to provide an 
adaptive control with higher precision, ensuring then low computational burden, suitable for high Dof robotic 
systems. 

The remainder of this paper is organized as follows: in Section 1, the dynamic model of each subsys- 
tem of n-DOF robot manipulators are described. In Section 3, the proposed TDVDC controller is designed, 
and compared to the conventional regressor-based virtual decomposition control is presented. Ine Section IV 
he proposed TDVDC is designed with stability analysis. Section V, a case study is performed on 7-DoF robot 
manipulator with TDVDC, and compared to the regressor-based VDC. The conclusions are summarized in 
Section VI. 


2. SYSTEM MODELING 
2.1. Desription of system 
The equation of motion of an n-DOF robot manipulators are described according to the Euler-Lagrange 
theory [3], as (1): 
7 — M(q)d + Clq, d) + G(q) (1) 


where M (q), C(q, d), G(q) and G(q, d) are respectively the manipulator's mass matrix, the Coriolis and cen- 
trifugal terms vector, the gravity terms vector, and the torque friction vector. In the VDC controller, the dynamic 
equations of the system can be expressed as link subsystems and joint subsystems, where Virtual Power Flow 
(VPF) [21] is used to characterize the coupling dynamic interactions among subsystems. 

The Figure | represents the virtual decomposition of serial robot manipulator to 7 links , where i = 
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1,...,.n, connected via mechanical joints. Each link has one driving cutting point according to the frame Bj, 
and one driven cutting point according to the frame B;. The it” joint has one driven cutting point according to 
the frame D; and one driving cutting point according to the frame T;. The dynamic equation of every subsystem 
is derived with respect to the local frame D; according to the Denavit Hartenberg formalism. 
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Figure 1. Virtual decomposition shematic of serial robot manipulator [21] 


2.2. Link dynamics 

The dynamic equation of the i*” rigid-link subsystem following its fixed frame can be expressed as 
(2) [21]: , 

met = Mg, V) + Cn, V Gs, Q) 

where P: F* denote the net force/moment vectors applied from the lower (i — p^ link to the i^^ link expressed 
in frame B;. P:V is the generalized linear/angular velocity, and Mp,, C p,, G p, represent respectively the 
inertial, Centrifugal/Coriolis, and gravitational terms, respectively. Using an iterative process computation, the 
vector of the total generalized force (forces/moments) acting on the i^^ rigid body can be computed as (3): 


Ba F mE Bn p* 
Bs-ip — Baa pe 4 Belts Bn p* 


Bip — Bi F* + B; Up, F* (3) 


where” F is the generalized force exerted by body i + 1 on body i, P: F is the generalized force exerted by 
bodyi + 1 on body i, P: V, is the velocity of body i. P:UT p, 44 İs the transformation matrix, defined as (4): 
Bi 

Hp, 03x3 (4) 


BirrT — 
U Biji = S Bi Bin BiR 
( Tax) B Bi+1 


i441 
where P: R B, ,,Iepresents the rotation matrix from frame B; to the frame B;,.1, 0323 is the null matrix, S is 
the skew-symmetric matrix operator performing the cross product between two vectors, and Pr, denotes a 
vector from the origin of frame D; to the frame B;1, expressed in frame B;. 
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2.3. Joint dynamics 
The dynamic equation of the 7 
(5) [21]: 


th joint subsystem expressed in its fixed frame is given by the following 


Tij = mide + kasign(di) (5) 
where Im; and ke; denote the moment of inertia, and the Coulomb friction coefficient of the i*” joint respec- 
tively. 

Finally the control torque of the global system can be expressed by (6): 


Ti = Tij t Til (6) 


Tij and T; denote the net torque and the control torque applied to the i” joint respectively. 7j; denotes the 
output torque of the it” joint toward the links, which can be computed in magnitude by the torque projected 
from the links, expressed in (7). 

mcg EF (7) 


where z = [0,0,0,0,0,1]” for revolute joint. d; the joint velocity vector. 


3. CONTROLLER DESIGN 

The main control objective in VDC approach is to track a required trajectory such that the joint track- 
ing error between the actual and required velocity converges asymptotically to zero in finite-time, with high 
accuracy even in presence of uncertainties and external disturbances. 
In order to design the VDC controller some vectors must be defined related to this approach. 


3.1. Required vectors 

The required velocity, is one of the concept related to VDC approach [21], wich can be expressed as 
(8): 

Gr = da + (da — q) (8) 

dq and qq denote respectively the desired joint velocity and the desired joint angle, A > 0 is a constant. 
The dynamic equations for the VDC controller design are based on the required joint velocity and the required 
linear/angular velocity vectors. 
The required linear/angular velocity of the link can be computed as (9): 


SHOE = xui PRU eas Ve (9) 


Adaptive control law of rigid-link subsystem 


32. Adaptive control law of link subsystem 
Referring to the dynamic link subsystem (2), and the required linear/angular velocity vector and its 
time-derivative, the required force/moment vectors are expressed as (10) [21]: 


d 
PF, = Mp7 V.) + OB,” V, Gs, (10) 


Where P: F,* the required net force/moment vectors of the subsystem links. 
PiV the vector of the generalized velocities (i.e., linear and angular components), wich can be expressed as 


(11): 
Boy oq perdu (11) 
Consider the linear parameterization form, the link subsystem (2) can be expressed as (12): 
Bi F,* = Yon (12) 


where the Y;; denotes the it” regressor matrix formed by the joint velocity, the linear/angular velocity and its 


time-derivative; and 6; denotes the i^^ parameters vector formed by the uncertainties parameter vector. 


Therefore, the control law of link subsystem is designed as (13): 


Bir = Yuĝu + Ku”'ey (13) 
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where Ky; is a diagonal matrix representing the gain of the feedback controller, and P:ey is a measure of the 
tracking accuracy defined by (14): 
Brey =?! V, By (14) 


Êi; is the estimate of the uncertainties parameter vector 0j; Finally, the control law of link subsystem, can be 
computed by an iterative process, as (15): 


DAE. = Yi Op, + Kin?" ey 


PF, = Yabu + Ku ey + Up,” (15) 


3.3. Adaptive control law of joint subsystem 
For the control law of joint subsystem defined as, the required net torque 7;;,. applied to the i*” joint, 
is based on the required joint velocity vectors, as (16) [21]: 


Tijr = Imi Gir + kei sign( dir) (16) 


According to the linear parameterization property, the required net torque 77. can be written in linear form as 
(17): 

Tijr = Yjiðji (17) 
where Yj; denotes the regressor matrix formed by the joint velocity and acceleration, and 0;; denotes the 
parameters vector formed by the physical dynamic parameters. 

Due to the difficulty in knowing the exact value of the physical parameters of the it” joint, they should be 
estimated. Then the estimation vector denoted by 654 is used, and the equation of control becomes (18): 


Tijr = Yjib; + Kjieq (18) 


where /€;; is a diagonal matrix representing the gain of the feedback controller, and e, is a tracking joint error 
defined by (19): 

€gi = dir — di (19) 
Finally, the total control torque is computed using the required output torque of the i*” joint toward the links, 
and the required control torque of the i^^ joint as (20): 


Ti = Tijv + Tilr (20) 


where 7;;,. denotes the control torque of the ?th joint, and T; the required output torque of the ¿th joint toward 
the links expressed with the required force/moment vectors as (21): 


Tilr = ge: (21) 


The control based VDC approach is to resolve equation (20), where the vectors of parameters estimation 61; 
and 0;; are used. The parameter adaptation function should be chosen to ensure system stability. 


3.4. Regressor-based VDC controller 


In the conventionnal Regressor-based VDC controller, the uncertainties parameter vectors Êj and 61; 
for joint and link subsystems, are updated using the projection function P defined as a differentiable scalar 
function [21]. According to the link subsystem, the uncertainties parameter vector is estimated as (22): 


ó;., = Psy; (t), Piy: Qiy (t), biy (t), t) (22) 
where 0j, denotes the yth element of Dis Si, (t) denotes the yth element of s;(t) defined as (23) [21]: 


silt) = Yj” (P V, -P: V) (23) 
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Piy > Qis a parameter update gain, and a; (t), biy (t) denote the lower and upper bounds of jiy. 
The projection function P is a differentiable scalar function defined by its time derivative which is governed by 
(24) and (25): 

P(t) = ps(t)k (24) 


with 


0 if P(t) < a(t) and s(t) « 0 
k= <4 0 ifP(t) > b(t) and s(t) >0 (25) 


1 otherwhise 


It is the same for i; of the required net torque. 

The use of the projection function for the uncertainties parameter vectors estimation, requires to com- 
pute the derivative of the regressor matrix in every sampling cycle. However, the regressor matrix derivation is 
not unique, though the process is standardized. Furtheremore, it presents high complexity, then an additionnal 
computational burden 


3.5. Proposed TDVDC controller 

As demonstrated, regressor-based VDC method presents inherent limits regarding the use of the pro- 
jection function to estimate the unknown parameter vectors 6; and 651 for link and joint subsystems respec- 
tively. To over come the issues with regressor-based VDC technique, a new control strategy combining the 
TDE and VDC approachs is proposed. The idea is refered to estimate the dynamic uncertainties and parameter 
vectors by the use of TDE. Refering to the dynamic link subsystem given in the (2), and the dynamique joint 
subsystem given in the (5), the dynamic uncertainties and unknown parameter vectors can be regrouped as (26) 
and (27): 
For the link subsystems 


Bi p* = Mg, — (PV) + Hu (26) 


For the Joint subsystems 
Tijr = Imidir + Hj (27) 


Hii and Hj; represents the dynamic uncertainties and unknown parameter vectors of the link and joint subsys- 
tems respectively, where (28) and (29): 


Hy = Cp, PV. + Gg, (28) 


Therefore, the control law subsystems, are given by (30): 


BiP * = m” V + Hi + KuPiev 
(30) 


Tijr = lidir + H ji + Kyieg 


m; and i; are a constant coefficients associated to M. B,, and Im; respectively. The determination of both con- 
stant coefficients mM; and 7; is discussed in [37], [38]. Hj; and Hj; represents respectively the estimate of Hj; 
and Hii. 
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In order to design the TDVDC controller and carry out its stability analysis, let consider the following 
assumptions: 
A1: The joint position and velocity are measured. 
A2: The parameter vectors H;; and Hj; 
their time derivatives åH 1, and 4H ji are globally Lipschitz functions. 
A3: The constant coefficients m; and i; are chosen assuming that: 
lI, — M(a)m | « 1 
=—1 
|In- I)i | «1 
According to the use of TDE [37], and if the assumption A2 is verified, we can estimate H); and H j;. 
Indeed the value of the function Hj; and H ;; are considered at the present time t, very close to that at time 
(t — T) in the past for a small time delay T' in (31). 


For the link subsystem . . 
Hy(t) = Hult — T) (31) 


therefore, using an iterative process, the estimate of the uncertainties parameter vector of the link substem 
H(t) can be computed as: 
Fi; (t) ES Ti (t = T)z - mP V(t = T) — KyP:ey(t ES T)- 
BiUp,. (t—- T) Ayisiy(t - T) + mi 4P^V(t- T)- 
Ka+” n ey(t — T) 
Fray) = (rixae (t — T)z — m4 PV (t - T)- 
Koa ey(t—T) -Pi Up, (t — T) Ayia) (t— T)+ 
Mio” tV (t — T) + Ki;,2P*?ey(t - T) 














Éi (t) ras (t — T)z — mn” V (t — T) — Ky, P^ey(t — T) 


For the joint substem, the estimates of the uncertainties parameter vector H ji(£) is given by (32) and (33): 


Hj(t) = Hj(t— T) (32) 
then : 
Hj(t) ~ Tijr(t =? T) = Kyiegilt — T) (33) 


where T is the estimation time delay. The accuracy estimation of Ñu; (t) and H,;(t) improves for a small T. 
In practice, the smallest estimation time delay T' is chosen to be the sampling period which means that the 
perfect parameters vector are identified every sampling period. Finally, the proposed control is obtained as 
(34)-(36): 
Tilt) = Tage (t) + Tar (t) (34) 
where : 
Tage (t) = fidis (t) + Hs (t) + Kyi (dir (t) — di(t)) (35) 


and 


Tar(t) = zT [m;P:V (t) + Ay (t) + Ky? ev (t) — 
Up. Mi V(t) + Braga (t) + Rigg ev (t)] (36) 


The closed-loop control system based on the proposed TDVDC technique is presented in Figure 2. 
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Figure 2. Block diagram of the proposed controller 


4. VIRTUAL STABILITY ANALYSIS 

According to the virtual work approach, the global stability of the system's VDC is proven through the 
virtual stability of each subsystem. Indeed, using the definition 2.17 and theorem 2.1 in [21], the global system 
is stable in the sense of Lyapunov, if each subsystem is proved to be virtually stable. It will be proven that all 
the decomposed subsystems of the studied system with their respective control equations are virtually stable, 
leading to the stability of the entire system. Generally, the stability analysis, in the sense of the Lyapunov 
approach, refers to define a positive candidate function and then to show that its variation is a decreasing 
function. Considering the Lyapunov candidate function for the entire robot as summation of two functions for 
the link (Vi;) and joint V;;)subsystems as (37): i = {0,...,n} 


V -M Wr Vj (37) 


4.1. Virtual stability of the ith link 
Let consider the Lyapunov candidate function for the i*” link as (38): 


1 1 i 
Vi = 5"ev Mp, "ey zu 5 Hu zi Hy ee) 


Then from [21] and the dynamic equation of the i^^ link given in (13), the first derivative along time of V;; can 
be given by (39): 


Vs = —Biey Ky, Bey + Biey’ (Bip, —Bi F) + 
(Hi; — Hu)(Piev — Hy) (39) 


where Bie, C BiP'ey = 0, since Cp; defined as skew-symmetric. 
According to the TDE use the V;; becomes (40): 


Vii = igs hee + Bi ey” (Bi Fe —Bi F,) + 
1 
AB (YT P:ey — 55 AHu) (40) 


where AHj(t) = Hu(t) — Hu(t — T), is the term due to the TDE error. 
Otherwise, as H);(t) is a Lipschitz function, then (41): 


[AH] € uT (41) 


du is the Lipschitz constant. 
To perform the VDC for each subsystem, the virtual power flows are introduced to characterize the dynamic 
interaction among the subsystems at its cutting points. Indeed, the virtual power flow is defined as the inner 
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product of the linear/angular velocity error vector and the force/moment error vector, with respect to the frame 
{A}, as (42): 
pa = tep (AF; —^ Fp) (42) 


Therefore from [21], (40)-(42), we obtain: 


. ru : 1 
Vi € —Piey KyPiey + pp, — pr, — 59 (43) 


where pg, and pr, represent the virtual power flows at the two cutting points of each link. 
As defined in [21], according to an open chaine structure, for pg, = 0 and pr,, = O the total virtual power 
flows is given by (44): 

5 (PB, — Pru) = 0 (44) 


Therefore the (43) becomes (45): 
X Vu < X (— Fey KuP'ey — E (45) 
: xg 2 


4.2. Virtual stability of the ith joint 
The positive Lyapunov candidate function related to the joint dynamics can be chosen according to 
the joint dynamic and its control law, as (46): 


1 1 E 
Vii = SInieq^ + 5 (Hji Hja)? (46) 
Then, its time derivative is (47): 
Vii = —egiIniégi — (Hji — Hyi) Hy (47) 


with the TDE use, and the dynamic equation of the it” joint given in (18), the Vi becomes (48): 
F * 1 
Vj; = — Kies + eqi(Tip — Tir) — gp AH (48) 
According to [21], (41), (48), and VPF definition, we obtain (49): 
. 1 
Vii < -Kjeh = 255i *DPB,; — PTj: (49) 
As described in the above section, using VFP the (49) becomes (50) 
; 1 
2 
»3 Vii < SC Ee = 35) (50) 


4.3. Stability of the global system 
The derivative of the global Lyapunov candidate function (37), is given as (51): 


V= s Ys 61) 


The V function is proved to be always decreasing based on the virtual power as the inner product of the 
linear angular velocity vector error and the force moment vector error presented in [21], and the choice of the 
parameter function adaptation, where (52): 


; NE 1 1 
V<- X (Pev Ku”tey F zii + KC aes + 355) (52) 
ij 
where ô; > 0 and à; > 0 are the Lipschitz constants. Since V < 0 where all gains are positive, the system is 


asymptotically stable in the sense of Lyapunov [21]. 
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5. CASE STUDY 
5.1. Simulation description 

To illustrate the effectiveness of the proposed control strategy, in this section a case study is performed 
for tracking trajectory of 7-DoF robotic manipulator using the proposed TDVDC. The simulation routine is 
conducted following the control architecture given in Figure 2, wich includes the desired trajectory given in 
joint space. The desired joint velocity and the desired joint acceleration are obtained from the derivation of the 
desired joint position.The equation of motion for each link and joint (i) (? = 1,...,7) subsystems is derived 
with respect to a local frame { B;} as shown in Figure 1. The mass, Coriolis, and the gravity termes of the link 
(i) can be described by (53)-(55): 


0 Mi —mdi 
Cp, = | mi 0 0 do t= Ay eel (54) 
misin(qi)g 
Gp, = | mi,cos(qi)g | , à — 1,...,7 (55) 
madicos(qi)g 


where the physical parameters of the using robot system are represented in Table 1. The numerical simulations 
are conducted for the proposed TDVDC controller and compared to the conventional regressor-based VDC 
in order to prove the effectiveness of the prposed approach. During the trajectory tracking, a disturbances 
was added to the torque input representing 5% of maximum value of the torque after t = 10s. In addition an 
uncertainty function U (t) was injected to the robot dynamic model to validate the effectiveness of the proposed 
control strategy in (56). 

U(qi,t) = qi * sin(t) + 0.5 x sin(500 x pi x t); (56) 


Table 1. Physical parameters 
l4 = 0.3m; l2 = 0.5m; l3 = l4 = 0.21m; l5 = 0.25 m; le = 0.5 m 
mı = 0.122Kg; m2 = 0.66K g; m3 = 0.08Kg; m4 = 0.175 Kg; ms = 0.251Kg; mg = 0.023 Kg 
key = key = keg = ke, = kes = keg = 0.5 N.m 
h = h = ÍI = I4 = Is = Ig = 0.0234 Kg.m? 























For the proposed TDVDC approach, the target robot is controlled following the closed-loop given in 
Figure 2. It concerns the use of TDE for the estimation terms defining the unknown and uncertainties parameter 
vectors of the robot. The required linear/angular velocity and its time derivative is computed using A constant. 
The constant coefficients mi; and i; are chosen according to the assumption A3. A suitable choose of these 
constants influence the stability and attenuation of measurement noise. These constants are conducted by the 
trial and error method. The time delay T is fixed as sampling time. The gain parameters of the feedback 
controller K; and K; for joint and link subsystems respectivelyare are fixed ensuring the stability condition. 
These parameters values must be adjusted in order to obtain the optimum performance. For the conventional 
regressor-based VDC approach, the parameters estimation is based on projection function presented in (22) 
which requires the derivation of the regressor matrix in every sampling time, as discussed previously. To 
accomplish the simulation routine, in addition to the gains feedback controller K}, K; and A, the parameters p, 
a, b are used for the projection function. 


52. Simulation results 

The obtained simulation results of the tracking trajectory and the traking errors for the proposed 
TDVDC and the conventional Regressor-based VDC strategies are shown in Figure 3, Figure 4, and Figure 
5 respectively 
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It is shown that the system response converges for the conventional regressor-based VDC and the 
proposed TDVDC approaches, with insensitivity to uncertainties and disturbances. However, for the proposed 
TDVDC the signal response present some delay at the begining before t = 0.2s. The obtained joint torque for 
the TDVDC and Regressor-based VDC are given respectively in Figure 6 and Figure 7. It should be noted that 
the feedback gains should carefully be selected to get an accurate tracking. 
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Figure 3. Tracking joint trajectory (a: for the joint q1, b: for the joint q2, c: for the joint qs, d: for the joint q4, 
e: for the joint qs, f: for the joint qe, g: for the joint q7) 
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Figure 4. Tracking error of TDVDC 


5.3. Discussion 

The goal of the proposed controller is to track the desired trajectory rather than convergence of physi- 
cally uncertain parameters. For the two numerical simulations conducted using the proposed TDVDC and the 
conventionnal regressor-based VDC, the robot can track the desired references but with relatively a tracking er- 
ror joint difference. In order to carry out an overall quantitative evaluation of the used controllers, the following 
evaluation criteria are adopted: Root mean square (RMS) of the tracking position error, defined by (57): 





Ti 
RMSļe:(t)] = y » f le;(t)2, dt (57) 


where e;(t) = qci — qai 
Computing time (CT) defined as a period from the start to the end moment of run simulation. 
Maximum torque (MT) defined as the largest torque input to the joint actuators. 

We can conclude that the proposed TDVDC presents the low tracking error compared to the con- 
ventionnal Regressor-based VDC. For the generated torque, the proposed controller present the low MT for 
the second joint q2, comparing to the regressor-based VDC for the same joint. In the otherhand, the CPU 
time with the proposed VDTDC approach is equals to 2.486009(s) and for the conventionnal VDC is equals 
to 4.721131(s), using MATLAB program. Therefore, in addition to the gain in terms of error the proposed 
VDTDC approach presents gain in terms of generated torque and also the computation time. 





error joint (rad) 














time(s) 


Figure 5. Tracking error of Regressor-based VDC 
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Figure 6. The obtained joint torques for the VDTDC controller 
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Figure 7. The obtained joint torques for the Regressor-based VDC controller 


Table 2. Comparison of RMS errors for regressor-based VDC and VDTDC. 








Joints Regressor-based VDC VDTDC 
qı 2.803404 0.13704 
q2 5.8692e704 4.0210—94 
qa 5.949994 1.201004 
qa 3.5726 04 3.569905 
q5 4.0829 04 4.6723-05 
q6 1.939494 1.3883- 05 
q7 5.2568 04 2.090404 





6. CONCLUSION 


To sum up, in this paper, a new aproach has been developed of nonlinear control of n-DoF robot 
manipulator. The proposed controller is a modified VDC introducing the TDE technique to estimate unknown 
and uncertainties vector parameters.The system modeling, and the VDTDC controller design were presented. 
Hence, the stability of overall system was proved based on the Lyapunov theory. Simulation results performed 
with 7-DoF robot manipulator for the proposed TDVDC controller, shows the capability of the controller for 
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tracking the predefined trajectory. Comparatively to the conventionnel regressor-based VDC, the TDVDC 
presents the lower tracking error, with the lower value of maximum torque, and finally lower computation 
time, wich proves the effectiveness of the propsal. As futur work, the proposed controller were be tested 
experimentally to another kind of application. 
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